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ABSTRACT 


The  performance  of  the  single-ended  multiangle  lidar  for  measurements  of  atmospheric 
extinction  and  backscattering  is  analyzed  in  the  presence  of  signal-induced  shot  noise.  An 
algorithm  for  determining  the  transmission  coefficient  in  the  presence  of  horizontal  inhomo¬ 
geneities  is  described.  The  analysis  indicates  that  real-time  measurements  of  atmospheric 
properties  with  the  multiangle  hdar  may  be  feasible  even  in  inhomogeneous  regions. 
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1.  Introduction 

Since  the  aerosol  constituents  of  the  atmosphere  can  be  highly  variable,  it  is  im¬ 
portant  to  have  the  means  to  determine  the  atmospheric  extinction  and  backscattering 
coefficients  in  real  time.  Remote  sensing  hdars  axe  the  most  promising  instruments  for 
such  measurements.^  Typically,  remote  sensing  lidars  require  ad  hoc  assumptions  about 
the  relationship  between  extinction  and  backscattering  in  order  to  determine  the  transmis¬ 
sion  coefficient  from  the  range-resolved  lidar  returns.  It  is  of  course  desirable  to  determine 
the  transmission  coefficient  from  lidar  data  without  assuming  any  such  relationships.  Sev¬ 
eral  methods  to  achieve  this  have  been  proposed:  the  Doppler  lidax,^’^  the  bipath  lidax,^’® 
and  the  multiangle  lidar.®“^° 

The  Doppler  lidar^’^  uses  differences  in  the  thermal  Doppler  broadening  of  radiation 
backscattered  by  molecules  and  aerosols  to  determine  their  respective  contributions  to  the 
overall  measured  lidar  return.  The  transmission  coefficient  is  then  calculated  by  compet¬ 
ing  the  measxired  molecular  contribution  to  that  predicted  by  the  standard  atmospheric 
models. 

The  bipath  method  employs  two  separate,  single-ended  lidar  systems  to  measure  both 
the  backscattering  and  extinction  coefficients  imambiguously.^’®  Hughes  zmd  Paulson®  have 
shown  that  the  backscattering  and  extinction  coefficients  obtained  from  a  single  lidar  using 
typical  a  priori  assumptions  (i.e,  Klett’s  reconstruction  algorithm^ were  significantly 
different  from  the  values  obtained  with  the  bipath  lidar.  However,  while  elegant,  the  bipath 
lidar  is  difficult  to  use  as  a  function  of  altitude. 

In  this  paper,  the  performance  of  a  single-ended  multiangle  hdar  for  real-time  mea¬ 
surements  of  atmospheric  transmission  and  beickscattering  is  discussed.  This  hdar  utilizes 
the  dependence  of  the  atmospheric  transmission  coefficient  on  the  zenith  angle,  to  ob- 
tmn  two  independent  equations  for  the  backscattering  and  transmission  coefficients  from 
measurements  along  two  paths  with  different  zenith  angles.  Therefore,  no  ad  hoc  assump- 
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tions  about  the  aerosol  extinction  and  backscattering  coefficients  are  required  in  order  to 
determine  these  coefficients  as  functions  of  zdtitude. 

Lidar  measurements  of  atmospheric  extinction  and  backscattering  as  a  function  of  the 
zenith  angle  have  been  carried  out  in  the  past  by  Spinhirne  ei  aV  Their  data  indicate  that 
at  the  time  of  measurements  the  atmosphere  was  horizontally  homogeneous  above  3  km, 
but  inhomogeneous  in  the  boundary  layer.  This  is  not  surprising  since  these  data  were 
collected  at  widely  spaced  intervals  of  the  zenith  angle  (which  was  varied  from  0°  to  about 
80°). 

Similar  measurements  were  attempted  by  Paulson,®  who  collected  lidar  returns  from 
optical  paths  with  elevation  angles  of  25  and  50°.  Even  though  the  atmosphere  was  not 
horizontally  homogeneous  for  such  large  angular  separations,  no  attempt  to  vary  the  ele¬ 
vation  angles  was  made. 

Successful  determinations  of  transmission  coefficient  from  the  lideir  returns  collected 
for  different  zenith  angles  were  reported  by  Sandford®  and  Shimizu  et  al.^°  These  results 
suggest  that  the  imcertainty  in  the  values  of  the  aerosol  extinction  and  backscattering 
obtained  by  Spinhirne  et  alJ  and  Paxilson*  could  be  considerably  reduced  by  utilizing  only 
the  returns  from  optical  paths  with  angular  separation  smaller  them  the  angvdair  size  of 
horizontal  inhomogeneities.  In  the  multiangle  lidar  discussed  here,  this  eingular  separation 
can  be  determined  in  real  time  by  comparing  returns  obtained  for  a  fixed  zenith  angle, 
but  different  a.zimuth  angles  of  the  optical  paths.  The  minimum  angular  separation  of  the 
paths,  based  on  noise  considerations,  is  discussed  in  Sec.  2.  An  algorithm  to  determine 
the  extinction  and  backscattering  coefficients  in  the  horizontally  inhomogeneous  region  is 
described  in  Sec.  3.  Results  of  calculations  using  this  algorithm  are  presented  in  Sec.  4 
and  conclusions  in  Sec.  5. 

2.  Multiangle  Lidar  Performance  in  a  Horizontally  Homogeneous  Atmosphere 

Consider  a  standend  lidar  pointed  along  a  LOS  with  zenith  angle  6.  The  lidar  equation 
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can  be  written  as^ 


E{z,9)  =  —  cos^  6  ^(z,6)  exp  — 


f  dz' kext{z' ,9)  , 
Jo 


where  E  is  the  average  energy  of  radiation  backscattered  from  the  range  R  =  z/  cos  9,  C 
is  the  calibration  constant,  and  P  and  k^xt  are  the  volume  backscattering  and  extinction 
coefficients,  respectively. 

Consider  two  lines-of-sight  with  zenith  angles  9  and  9  +  A0,  respectively.  Assuming 
that  the  atmosphere  is  horizontally  homogeneous  over  the  angular  separation  A0,  the  two- 
way  transmission  coefficient  as  a  function  of  range  can  be  obtained  from  the  lidar  signals 
obtained  for  zenith  angles  9  and  9  +  A9 

T{R  =  z /  cos  9)  =  exp - f  dz' kextiz\0)  =  [^(z)]'^  (2) 

L  cos  9  Jq 


where 


cos(9  -f  A9) 
cos(9  +  A9)  —  cos  9  ’ 


^  E{zl  cos  9)  cos^(g  -I-  A9) 

E{z/  cos{9  +  A9))  cos^  9 

Since  Q{z)  depends  only  on  the  ratio  of  the  received  signals,  the  transmission  coefficient 
can  be  determined  without  absolute  calibration  of  the  lidar.  This  is  also  the  case  for  the 
extinction  coefficient  which  is  given  by 

Kxt{z)  =  -\i  ^  1“  Q{^)  ■  (5) 

However,  absolute  calibration  of  the  lidar  is  still  required  in  order  to  determine  the  voliime 
backscattering  coefficient  ^  from  the  lidar  equation  (1). 

Since  7  increases  with  decreasing  angular  separation  between  the  two  optical  paths, 
there  must  be  a  minimum  value  of  A9,  or  maximum  value  of  7,  in  order  to  obtain  reliable 
results  in  the  presence  of  noise.  The  performance  degradation  of  the  multiangle  lidar  due 
to  signal-induced  shot  noise  is  discussed  below. 


FR-0215 


EOS 


Let  Nk{z,  9)  be  the  number  of  photocounts  detected  by  the  lidar  receiver  from  altitude 
z  along  the  optical  path  with  zenith  angle  9  for  the  A:th  pulse.  If  turbulence  effects  axe 
neglected,  the  probability  of  detecting  such  a  number  of  photocounts  is  governed  by  Poisson 
statistics,  i.e.,  the  probability  density  function  is 

-  (n.+N,)  /0\ 

where  Ng  and  Ni,  are  the  expected  numbers  of  signal  and  background  photocounts,  respec¬ 
tively.  The  expected  value  and  variance  of  the  detected  photocounts  are: 


<Nk  >^N,  +  Nb, 


(7) 


ah  =  <  Nl  >  -  <  Nk  >^  =  Ns  +  Nb. 


(8) 


In  order  to  reduce  the  variance  in  the  measurements,  consider  a  sum  of  random  variables 
Nk: 

1  A 

(9) 


fc=l 


where  K  is  the  nimiber  of  laser  pulses  used  in  averaging.  Then,  for  iiC  — »  oo,  the  central 
limit  theorem  gives  the  probability  density  function  for  n-}^ 

1^/2  r  i\T  ^r.^2^ 


p(n)  = 


K 


[2ii^{Ns  +  Ni,)\ 


exp 


-K 


jn-Ns-Nb) 

2{Ns  +  Nb) 


(10) 


The  corresponding  expectation  value  and  variance  are: 


<n>  =  Ns  +  Nb,  (11) 

al^{N,  +  Nb)/K.  (12) 


The  purpose  of  the  measurements  is  to  determine  the  two-way  transmission  coeflficient 
from  the  groimd  to  altitude  z  along  the  optical  path  with  zenith  angle  9.  Let  ni  and  n2  be 
the  average  numbers  of  photocounts  from  adtitude  z  for  the  opticzd  paths  with  zenith  angles 
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6  +  Ad  and  6,  respectively.  Then,  assuming  that  horizontal  gradients  in  baekscattering 
and  extinction  coefficients  can  be  neglected,  the  transmission  coefficient  is 


where 


and 


T{z,d)=A'^ 


A  = 


n2  -  Nb2 
cos^(d) 


7  = 


cos^(d  +  Ad) 
cos(d  +  Ad) 


(13) 


(14) 


(15) 


cos(d)  —  cos(d  +  Ad) 

Since  n\  and  n2  axe  independent  random  variables,  the  expectation  value  and  vciriance 
of  T  can  be  written  as: 


<T>  =  A'' Nbi ,  K, )  I.y{N,2  ,Nb2,K2), 


(16) 


and 


4  =  A^-'[h,{N,„Ni,uKi)I-2.,(N.2,Nu,K2)-JUN,uNu,Kt)ll^(N.2,NK,K2)], 


(17) 


where 

l4N„Nt,K)  = 


K 


L2ir(Ar,  +  Nt) 


^OO 

/  dn{n-  Nb)°  exp 


-K 


{n-Ns-NbY 

2{N,  +  Nb)  J 


(18) 


A  small  positive  vzJue  of  €  (e.g.,  e  =  0.01)  is  introduced  in  order  to  avoid  artificial  sin- 
guleuities  for  negative  values  of  a  if  n  =  iVj.  Cleeirly,  in  the  absence  of  signal  no  useful 
information  about  the  transmission  coefficient  can  be  obtmned  from  the  measurements. 


It  should  be  noted  that  the  range  of  integration  in  Elq.  (18)  is  such  that  the  veJues 
of  n  <  Nb  +  eNa  are  excluded.  This  is  a  very  good  approximation  since  the  probability 
of  n  in  this  range  is  less  than  0.5%  for  KNal{Na  +  Nb)  >  5.  As  will  be  seen  later,  the 
condition  KNg/{Na  +  Nb)  ^  5  must  be  satisfied  in  order  to  obtain  reliable  values  of  the 
transmission  coefficient. 
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In  general,  the  integrals  in  Eq.(17)  must  be  evaluated  numeric2illy.  However,  in  the 
limit  K  — ♦  oo,  these  integrals  can  be  estimated  analytically.  Eq.  (18)  may  be  rewritten 

N°  f°°  /  j 

Ia{N,,Ni,K)  =  (1  +  t)  e-^  (19) 

(20) 

For  ^  1  we  have 


where 

.2  _ 

2(N,  +  Nty 


On  =  0(0:  —  1)  •  •  •  (a  —  n  +  l)/n! ,  n>l.  (22) 

Therefore,  the  relative  error  in  the  value  of  the  transmission  coeflBcient  obtained  with  the 
multiangle  lideir  is 

O'r  _  2  ^9l  +  N(,i  Na2  +  Ni,2 

<T>^~'’  t  KiN]j  KiN],  . 


I  7"(T-l)(37-5)(iV.i+JVM)^ 


7%  +  l)(3T  +  5)(Jf.;  +  iVt2)" 

2kin;, 


(23) 


+ 


37<(JV„+iVu)(iV.2  +  JV„) 
2KlK^Nj,Nj^ 


+  0(S-'‘). 


Expression  (23)  is  useful  for  estimating  the  number  of  laser  pulses  Ki  and  K2  required 
to  obtain  the  desired  accuracy  in  transmission  measurements,  given  the  signal  and  back¬ 
ground  levels  along  the  two  optical  paths.  It  can  edso  be  used  to  determine  the  maximum 
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allowed  value  of  7,  i.e.,  the  minimum  angular  separation  between  the  optical  paths,  if  other 
pzurzuneters  are  fixed. 

The  ch2iracteristic  length  scale  of  horizontal  inhomogeneity  can  be  determined  by  the 
multiangle  lidar  in  real  time  by  collecting  returns  for  a  fixed  zenith  angle  and  different 
azimuth  angle.  For  a  horizontally  homogeneous  atmosphere  the  lidar  return  should  be 
independent  of  the  azimuth  angle  (f>.  In  practice,  this  can  be  tested  by  comparing  the 
range-resolved  lidar  returns  from  optical  paths  separated  by  A<f>,  taking  into  account  the 
uncertmnties  of  individual  returns.  If  the  atmosphere  is  not  homogeneous,  the  separation 
in  eizimuth  angle  could  be  decreased,  until  the  locally  homogeneous  region  has  been  found. 
If  the  corresponding  separation  in  the  zenith  angle  is  less  than  the  minimum  sepairation 
allowed  by  the  noise  considerations,  then  it  may  become  necessary  to  use  a  different  ap¬ 
proach.  An  algorithm  designed  to  perform  in  the  presence  of  inhomogeneities  is  described 
in  the  following  section. 

3.  An  Algorithm  to  Determine  Extinction  and  Backscattering  in  a  Horizontally 
Inhomogeneous  Atmosphere 

In  the  derivation  of  the  transmission  coefficient  from  the  hdar  returns  for  two  optical 
paths  with  zenith  angles  0  and  0  -f-  A0  it  has  been  assumed  that  the  backscattering  and 
extinction  coefficients  are  functions  of  altitude  only.  Consider  now  a  situation  when  this 
assumption  is  not  valid. 

Prom  Eq.  (1),  for  a  fixed  zizimuth  angle,  the  rzinge-corrected  lidar  return  from  the 
optical  path  with  the  zenith  angle  0  is 

F{z,0)  =  CP(z,0)  exp 

=  C^iz,0)T{z,0). 

In  general,  the  backscattering  and  extinction  coefficient  can  be  written 

P{z,  0)  =  /3o(^)  [1  +  fp{z,  0)] ,  (25) 


COS0 


f  dz' 
Jo 


ktxtiz',0) 


(24) 
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and 

h,rt{z,e)=ko{z)[l+fk{z,e)],  (26) 

where  fp  and  fk  are  nonnegative  and  vanish  in  the  absence  of  horizontal  inhomogeneity. 
If  mecisurement  of  backscattering  and  extinction  is  desired  for  some  zdtitude  z  and  zenith 
angle  then  it  is  clear  that  the  performance  of  the  multiaingle  lidar  could  be  considerably 
improved  if  6j„,  the  zenith  angle  of  the  second  optical  path,  could  be  chosen  to  minimize 
the  effect  of  horizontal  inhomogeneities,  i.e.,  if 


(27a) 


and 


fk(^,0m)  =  h(‘,0n) 


(27b) 


Then,  the  transmission  coefficient  can  be  determined  from 


T(z,dn)  = 


F{z,dn) 

[Fiz,0m) 


Tmn 


(28) 


where 


Ifmn  — 


COs{dm) 


(29) 


COs{0m)  -  COS(0„)  ■ 

Eqs.  (27)  cannot  be  solved  simultaneously  without  additional  aissumptions.  Nevertheless, 
it  is  still  possible  to  choose  the  optimal  2mgular  separation  between  the  two  optical  paths; 
an  algorithm  to  do  so  is  described  below. 


1.  Lidar  returns  are  obtained  at  zenith  angles  dn  =  9o+nA9^  where  the  total  range  of  angles 
covered  should  be  greater  than  the  angular  size  of  inhomogeneity  under  consideration. 
Data  collection  should  include  averaging  over  a  number  of  pulses  eis  well  eis  backgroimd 
measurements  to  determine  the  expected  VEdue  Euid  standard  deviation  of  the  bEickscattered 
signal.  Prom  these  data  70,  the  mEiximum  value  of  7  consistent  with  the  desired  acciiracy, 
must  be  determined  using  Eq.  (23). 
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2.  In  the  absence  of  horizontal  inhomogeneity 


^  2  tan0 

—  In  F(z,d)  = - T- 

dd  ^  ^  COS0 


Therefore,  for  a  fixed  altitude, 


f  dz' kexi{z') 
Jo 


(30) 


0(6.)  =  ^ 


tan  9„ 


(31) 


if  the  atmosphere  is  horizontally  homogeneous.  The  homogeneous  part  of  the  lidar  return 
can  be  determined  by  finding  the  values  of  i  and  j  such  that 


|G(0i)  -  G(0j)l  =  min,  and  |7._,|  <  70 


In  order  to  avoid  accidental  minima,  the  value  of  i  can  be  chosen  from 


(32) 


|G(«i)  -  G(9i-,)1  +  |G(9i)  -  G(«i+,)|  =  min, 


(33) 


in  order  to  ensure  that  is  in  the  homogeneous  region.  Then  the  appropriate  value  of  j 
is  chosen  using  Eq.  (32).  The  homogeneous  part  of  the  lidar  return  is  given  by 


ft(z,«)  =  C,8»(z)[T»(2)|'/'”'  «  C0<,{z)[%(z)YI 


costf 


(34) 


where 


and 


Th(z)  =  [F(zA)m=^,«i)V» 


CMz)  =  F(z.9i)(r4(^)) 


— 1/  cos  Oi 


3.  The  next  step  requires  calculation  of  the  deviation  of  the  mezisured  lidar  return  from 
its  homogeneou  part  at  {z,0n)  by  determining 

AF(z,0„)  =  ln[F(z,5„)/Fft(z,0„)] 

(36) 


ln(l  +*(z,9„)l-  ^  [  dz'h{z')Mz',e„). 
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Two  limiting  cases  should  be  considered.  In  the  first  case,  the  meaisured  lidax  return  is 
dominated  by  increased  backscattering  and  AF  >  0.  In  this  case,  condition  (27a)  can  be 
satisfied  if 

/^F{z,e,n)  =  ^F{z,en).  (36) 

In  the  second  case,  the  measured  return  is  dominated  by  decreased  transmission  and  AF  < 
0.  Therefore,  to  satisfy  condition  (27b)  one  should  have 


cos  Ojn  AF{z,  drn)  =  COS  0„  AF(z,  6„). 


(37) 


Depending  on  the  sign  of  AF,  Eq.  (36)  or  (37)  is  used  in  the  algorithm  described  here  to 
determine  the  optimal  angular  separation  between  the  optical  paths  such  that  j^n  ^  7o- 
Then,  the  transmission  coefficient  for  each  angle  is  found  using  Eqs.  (28)  and  (29),  subject 
to  the  constraint  T(z,0„)  <  Tq(zY^^°^^’' .  If  AF  «  0,  then  the  transmission  coefficient 
along  that  LOS  is  the  same  as  in  the  homogeneous  case.  If  AF  has  an  extremum  for  some 
value  of  On,  so  that  Eq.  (36)  or  (37)  cannot  be  satisfied,  then 


T{z,en) 


F{z,  6n+i) 


1  7e« 


(38) 


where 


7e*  ~ 


cos(0„4.i)  COs{dn-j) 


<70- 


(39) 


cos(0„)  [cos(0„_j)  -  cos(0n+i)] 

The  values  i  and  j  in  Eqs.  (38)  and  (39)  can  be  found,  for  example,  as  follows.  For  j  =  1 
find  the  value  of  i  such  that  either  Eq.  (36)  or  (37)  is  satisfied,  subject  to  the  constraint 
given  by  Elq.  (39).  If  this  value  of  i  is  greater  than  some  specified  maximmn  value  (10 
was  used  in  the  calculations  described  here)  then  the  process  is  repeated  with  increasing 
values  of  j  imtil  all  of  the  above  conditions  are  satisfied. 


In  the  presence  of  a  horizontal  inhomogeneity  of  finite  extent  in  altitude,  say  from 
zq  to  zi,  the  algorithm  is  expected  to  perform  best  for  z  <  zo  and  z  >  zi.  The  worst 
performance  may  be  expected  if  the  change  in  lideir  return  due  to  increased  backscattering 
is  nearly  compensated  by  decreased  transmission. 
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Another  source  of  error  may  be  due  to  the  finite  angular  spacing  between  the  opticed 
paths  used.  The  return  at  the  optimal  value  of  the  zenith  angle  may  not  be  available.  This 
may  be  especially  important  for  inhomogeneities  of  small  angular  size. 

4.  Model  Calculations 

Algorithm  performance  was  tested  on  two  simple  models  of  horizontal  inhomogeneities: 

(A)  A  Gaussian  model  of  inhomogeneity 

f  exp  [-(«  -  ,  20  <  2  < 

m  =  \  (40) 

y  0 ,  otherwise 

(B)  A  double  Gaussian  model  of  asymmetric  inhomogeneity  with  structure  on  smaller 
angular  scale 

(  exp  [-(0  -  6cf/26l,]  +  0.5  exp  [-{6  -6c  +  40u,)^/20y  ,  zq  <  z  <  zi 

m  =  {  (41) 

y  0 ,  otherwise 

In  order  to  reduce  the  numbers  of  arbitrary  parameters,  it  is  assumed  here  that  ffi{z,  6)  = 
Ap  f{d)  zmd  fk{z,  6)  =  Ak  f{S).  It  is  important  to  note  that  this  simplification  does  not 
affect  the  performance  of  the  algorithm,  since  it  does  not  rely  on  any  prior  information 
about  the  structure  of  the  inhomogeneous  region,  which  is  usually  imavailable  during 
measmements. 

For  the  horizontal  inhomogeneity  extending  in  altitude  from  zo  to  zi,  let  To{z)  = 
exp[—  dz'  ko{z')]  and  Ti(z)  =  exp[—  dz'  Then,  Eq.  (24)  may  be  written  as 

F{z,e)  =  )9(z)[l  +  .  (42) 

The  following  cases  were  considered:  (1)  below  inhomogeneity;  (2)  at  the  lower  edge  of 
inhomogeneity;  (3)  and  (4)  inside  inhomogeneity;  and  (5)  above  inhomogeneity.  Suppose 
that  a  suflicient  number  of  pulses  is  used  for  signal  averaging  so  that  the  relative  error  in 
transmission  coefficient  from  Eq.  (23)  is  less  than  1%  for  70  =  15.  Then,  the  expected 
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lidajr  signal  for  the  two  models  is  shown  in  Figs.  1  and  2.  It  is  seen  that  at  the  lower 
edge  of  inhomogeneity  (case  2)  the  signal  is  dominated  by  the  increase  in  backscattering; 
in  case  4  the  increase  in  backscattering  is  almost  completely  compensated  by  decrease  in 
transmission  coefficient;  while  in  case  5  the  signal  is  dominated  by  reduced  transmission. 
The  performance  of  the  algorithm  for  these  lidar  signeds  is  illustrated  in  Figs.  3-4  and  5-6, 
where,  irrespective  of  the  sampling  interval  used,  the  calculated  transmission  coefficient  is 
shown  at  intervals  of  5°  to  avoid  symbol  overcrowding.  As  expected,  the  performance  is  the 
worst  for  case  4  where  the  lidax  signal  closely  resembles  that  of  homogeneous  atmosphere. 
The  performance  is  very  good  at  both  the  lower  and  upper  edges  of  the  inhomogeneity, 
indicating  that  this  method  can  provide  accurate  data  on  transmission  through  inhomo¬ 
geneities  such  as  clouds,  provided  they  are  not  optically  thick.  This  conclusion  appears  to 
be  independent  of  the  shape  and  nature  of  inhomogeneity;  the  results  are  equally  good  for 
both  models.  The  accuracy  does  depend  on  the  sampling  interval  Ad  used  during  mea¬ 
surements.  Not  surprisingly,  comparison  between  Figs.  3-4  and  5-6  indicates  that  more 
complex  structures  of  inhomogeneity  require  finer  sampling  intervals.  For  a  given  sampling 
interval,  the  accureicy  does  not  depend  strongly  on  the  value  of  70,  provided  it  is  chosen 
to  allow  the  use  of  neighboring  optical  paths  to  determine  the  transmission  coefficient. 

5.  Conclusions 

The  performance  of  the  single-ended  multiangle  lidax  for  real-time  measurements  of 
atmospheric  properties  was  amedyzed  in  the  presence  of  signal-induced  shot  noise  and  hori¬ 
zontal  inhomogeneities  in  atmospheric  extinction  and  backscattering.  Noise  considerations 
lead  to  the  minimum  angular  separation  between  the  optical  paths  (based  on  signal  and 
background  levels,  and  on  number  of  pulses  used  for  signal  averaging)  needed  to  deter¬ 
mine  the  transmission  coefficient  with  a  desired  accuracy  from  the  range-resolved  lidar 
returns.  A  simple  algorithm  to  determine  the  transmission  coefficient  to  a  given  altitude 
in  a  horizontally  inhomogeneous  region  was  described.  This  algorithm  makes  no  ad  hoc 
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assumptions  about  the  relationship  between  the  backscattering  ajid  extinction  coefficients, 
and  does  not  require  a  priori  knowledge  about  the  structure  of  inhomogeneities.  While  this 
algorithm  does  not  perform  very  well  inside  a  horizontzJ  inhomogeneity,  it  can  determine 
the  overall  transmission  of  optically  thin  inhomogeneities  with  high  accuracy. 

The  author  wishes  to  thank  M.  Elbaum  and  M.  Greenebaum  for  helpful  discussions. 
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Fig.  1.  Expected  lidar  return  for  a  model  A  of  horizontal  inhomogeneity  (Eq.  (23)): 

(1)  below  inhomogemeity  (To  =  0.6,  Ti  =  1.0,  Ap  =  0.); 

(2)  at  the  lower  edge  of  inhomogeneity  (To  =  0.5,  Ti  =  1.0,  Ap  =  1.); 

(3)  inside  inhomogeneity  (To  =  0.4,  Ti  =  0.8,  Ap  =  1.); 

(4)  inside  inhomogeneity  (To  =  0.3,  Tj  =  0.6,  Ap  =  1.); 

(5)  above  inhomogeneity  (To  =  0.2,  Ti  =  0.6,  Ap  =  0.). 
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Fig.  2.  Expected  lidar  return  for  a  model  B  of  horizontal  inhomogeneity  (Eq.  (24)); 
different  cases  are  defined  in  Fig.  1. 
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Fig.  3.  Algorithm  performance  for  model  A  of  horizontal  inhomogeneity  with  =  5®; 
different  cases  axe  defined  in  Fig.  1. 


FR-0215 


16 


ICOS 


Fig.  4.  Algorithm  performance  for  model  A  of  horizontal  inhomogeneity  with  =  1"; 
different  cases  are  defined  in  Fig.  1. 
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Fig.  5.  Algorithm  performance  for  model  B  of  horizontal  inhomogeneity  with  A9  =  5"; 
different  cases  are  defined  in  Fig.  1. 
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Fig.  6.  Algorithm  performance  for  model  B  of  horizonted  inhomogeneity  with  A0  =  1® 
different  cases  are  defined  in  Fig.  1. 


FR-0215 


19 


EOS 


References 

1.  R.  M.  Measures,  Laser  Remote  Sensing,  Wiley,  New  York,  1984. 

2.  S.  T.  Shipley,  D.  H.  Tracey,  E.  W.  Eloranta,  J.  T.  Trauger,  J.  T.  Sroga,  F.  L.  Roesler, 
and  J.  A.  Weinman,  “High  Spectral  Resolution  Lidar  to  Measure  Optical  Scattering 
Properties  of  Atmospheric  Aerosols.  1:  Theory  and  Instrumentation”,  Appl.  Opt., 
22,  3716-3724  (1983). 

3.  J.  T.  Sroga,  E.  W.  Eloranta,  S.  T.  Shipley,  F.  L.  Roesler,  and  P.  J.  Tryon,  “High 
Spectral  Resolution  Lidar  to  Measure  Opticzd  Scattering  Properties  of  Atmospheric 
Aerosols.  2:  Calibration  and  Data  Analysis”,  Appl.  Opt.,  22,  3725-3732  (1983). 

4.  G.  J.  Kunz,  “Bipath  Method  as  a  Way  to  Measure  the  Spatial  Backscatter  and  Ex¬ 
tinction  Coefficients  with  Lidar”,  Appl.  Opt.,  26,  794-795  (1987). 

5.  H.  G.  Hughes  and  M.  R.  Paulson,  “Double-Ended  Lidair  Technique  for  Aerosol  Stud¬ 
ies”,  Appl.  Opt.,  27,  2273-2278  (1988). 

6.  J.  D.  Spinhime,  J,  A.  Reagan,  and  B.  M.  Herman,  “Vertical  Distribution  of  Aerosol 
Extinction  Cross  Section  and  Inference  of  Aerosol  Imaginary  Index  in  the  Troposphere 
by  Lidar  Technique”,  J,  Appl.  Meteor.,  19,  426-438  (1980). 

7.  M.  R.  Paulson,  “Atmospheric  Horizontal- Inhomogeneity  Effects  on  the  Optical  Depths 
Determined  by  the  Double- Evaluation- Angle  Lidar  Technique”,  Report  TD  1600, 
DTIC  AD-a212  836,  NOSC,  1989 

8.  M.  C.  W.  Sandford,  “Laser  Scatter  Measurements  in  the  Mesosphere  and  Above”,  J. 
Atmos.  Terr.  Phys.,  29,  1657-1662  (1967). 

9.  H.  Shimizu,  Y.  Sasano,  H.  Nakane,  N.  Sugimoto,  I,  Matsui,  and  N.  TaJteuchi,  “Large 
Scale  Laser  Radsur  for  Measuring  Aerosol  Distribution  Over  a  Wide  Area”,  Appl.  Opt., 
24,  617-626  (1985). 

10.  J.  D.  Klett,  “Stable  Analytical  Inversion  Solution  for  Processing  Lidar  Returns”,  Appl. 
Opt.,  20,  211  (1981). 

11.  J.  D.  Klett,  “Lidar  Inversion  with  Variable  Backscatter  Extinction  Ratios”,  Appl. 
Opt.,  24, 1638  (1985). 

12.  W.  Feller,  An  Introduction  to  Probability  Theory  and  Its  Applications,  Wiley,  New 
York,  1970. 


FR-0215 


20 


